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In  a  recent  paper  White  and  McNally  (1987,  here¬ 
after  WM)  developed  and  discussed  a  particular  so¬ 
lution  to  low  frequency,  wind-forced  motion  in  a  strat¬ 
ified  ocean  of  constant  Coriolis  parameter/,..  There  is 
a  mixed  layer  of  depth  H.  where  the  horizontal  veloc¬ 
ities  U,„  I'.and  pressure,  P„.  are  independent  of depth; 
the  turbulent  stress  and  vertical  velocity  are  a  linear 
function  of  depth.  Below  is  a  stratified  deep  layer  of 
constant  Brunt-Vaisalii  frequency  N.  where  the  tur¬ 
bulent  stress  vanishes.  The  wind  forcing  is  a  sinusoidal 
stress  r  =  [0.  r1],  where  r  ‘  is  the  north-south  com¬ 
ponent  of  frequency  w  and  wavenumber  k  which  travels 
eastward  at  constant  phase  speed,  C\  =  w/k.  White 
and  McNally  present  a  solution  in  which  they  assumed 
ab  initio  that  pressure  gradients  in  the  mixed  layer 
vanish  and  then  force  the  lower  stratified  layer  with  a 
time-dependent  Ekman  suction  (which  is  independent 
of  mixed  layer  pressure  gradients).  Then  they  compute 
a  pressure  gradient  in  the  lower  layer,  and  assume  that 
this  lower  layer  pressure  gradient  also  acts  in  the  mixed 
layer,  producing  a  geostrophically  balanced  flow.  Such 
a  method  of  solution  is  incorrect,  as  it  is  self-contra¬ 
dictory  ( first  there  is  no  pressure  gradient  and  then  one 
is  arbitrarily  inserted  later)  and  it  is  incorrect  to  sup¬ 
press  any  pressure  gradient  in  fluid  motion  ab  initio, 
as  pressure  is  simply  a  reaction  to  conservation  of  vor- 
ticity  in  a  fluid  that  is  incompressible. 

In  this  note  we  presenWhe  correct  solution  to  the 
problem  WM  set  out  to  solve.  Qualitatively,  none  of 
the  results  presented  by  WM  need  be  changed.  Quan¬ 
titatively.  significant  differences  occur,  especially  as  w 
approaches  the  inertial  frequency  /,. 

With  pressure  terms  included,  the  expression  for 
mixed  layer  velocities  are  [Eq.  (3.1)  in  WM],  in  the 
usual  fluid  mechanical  notation. 
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In  the  above.  H’„  is  the  vertical  velocity  at  c  =  -//. 
and  in  its  derivation  we  have  used  the  condition  that 
IT  =  0.  at  r  =  0.  the  surface. 

The  vorticity  equation  in  the  stratified  layer,  as  ex¬ 
pressed  in  terms  of  pressure  is  [Eq.  (2.2)  in  WM]. 

£{(aT+£10  +  A'!v”!'p}-0'  ,4> 

and  IT  is  computed  from 

P.'H'  =  tt;  P.  (5) 

!\ '  dtdz 

The  horizontal  velocities  in  the  stratified  layer  have 
expressions  identical  to  those  in  the  mixed  layer,  except 
7=0.  The  critical  difference  between  ( 3 )  and  (3.1)  in 
WM  is  the  expression  in  the  last  term  on  the  right- 
hand  side  of  Eq.  (3). 

The  wind  stress  is 

r  =  [0,  r  ‘  Re(e,A '“'“')]  (6) 

and  following  WM.  we  are  interested  in  the  solution 
of  the  form 


P  =  R  e(P(z)e" 
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In  the  stratified  layer,  the  complex  valued  expression 
for  pressure  amplitude  that  satisfies  (4)  is 

P  =  Paeni!+,,),  (8) 

with 

,  n2  =  k2N7/(f02  -  w2 )  for  a  <f0.  (9) 

The  above  expression  matches  the  mixed  layer  pressure 
at  z  -  —H  and  pressure  perturbation  vanishes  as  z  -*■ 
-oo.  Substituting  (8)  into  (5)  and  evaluating  it  at  z 
-H,  where  W (z  =  -H)  =  Wa,  yields 

iun  - 

-jfi -Po=PoWo.  (10) 

The  second  relationship  between  Pa  and  H  '„  is  obtained 
from  (3),  (6)  and  (7);  it  is 

,  rl  2  \  tkf0T'  I  ?  r,  /  i  i\ 

Poifo  ~  w 2)  —  =  — — - u»k2P„.  (11) 


Eliminating  Pa  between  (10)  and  (11)  yields  the 
expressions  for  the  mixed  layer  velocity  response 
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Fig.  1.  Normalized  velocity  response  p„faIIL  „/t'  (a)  and 
-ip„/,//F„/r'  (b).  The  dashed  curve  is  the  solution  presented  by 
WM  and  the  solid  curve  is  the  correct  solution. 


Fig.  2.  The  difference  between  the  correct  solutions  for  curves  (a) 
in  Fig.  I  (solid)  and  curves  (b)  in  Fig.  1  (dashed). 


with 

ao2  =  kHN(f2  -  o>2),/2,  (15) 

Ro2  =  Poifo2  -«32  +no2).  (16) 


These  solutions  differ  from  those  presented  by  WM  in 
the  expression  for  R2.  If  the  last  term  on  the  right- 
hand  side  of  ( 3 )  is  a  priori  set  equal  to  zero  incorrectly, 
then  Ru2  =  pd  f,,2  ~  o>2).  WM’s  point  was  that  terms 
proportional  to  ft,2  are  important  in  the  expressions 
and  the  effect  of  the  pressure  gradient  in  the  mixed 
layer  as  w  -►  0. 1  f0  results  in  a  downwind  component 
of  velocity  in  Eq.  ( 14)  by  the  term  proportional  to  tt2. 

Figure  I  shows  the  correct  response  function  for 
PofoHUJH  and  —ip0f>HV0/ ry  as  a  function  of  a >//, 
with  HN/CX  =  tt/6  (the  same  value  used  by  WM), 
together  with  WM  solution.  Figure  2  shows  the  differ¬ 
ence  between  the  two  sets  of  curves.  Qualitatively,  the 
so'  ztions  do  not  differ  because  the  large  inertial  re¬ 
sponse  occurs  near  the  inertial  frequency  in  the  pa¬ 
rameter  range  used  by  WM,  regardless  of  the  pressure 
correction.  Furthermore,  given  a  W  at  z  =  —  //,  the 
physics  in  the  region  below  the  mixed  layer  is  treated 
correctly,  which  also  represents  a  portion  of  a  quan¬ 
titatively  correct  solution  and  contributes  to  the  qual¬ 
itative  agreement.  But  we  think  it  valuable  to  publish 
the  hydrodynamically  correct  response  function,  and 
to  give  the  correct  physical  interpretation. 

Acknowledgments.  This  research  was  sponsored  by  1—1  • 

NORDA’s  Tactical  Oceanography  Program  (element--  " 
62435N).  ‘ 


REFERENCE 

While,  W.  B..  and  G.  McNally.  1987:  Evanescent  pressure  gradient  ~ 

response  in  the  upper  ocean  to  subinertial  wind  stress  forcing  Codes 
of  finite  wavelength.  J  Phys.  Oceanogr..  17,  1032-1047.  d/or 


Dist 


Special 


m *l 


